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ABSTRACT 

We present a detailed analysis of the detached eclipsing binary CD Tau. A large 
variety of observational data, in form of IR photometry, CORAVEL radial velocity 
observations and high-resolution spectra, are combined with the published light curves 
to derive accurate absolute dimensions and effective temperature of the components, 
as well as the metal abundance of the system. We obtain: Ma = 1.442(16) M 0) 
R A = 1.798(17) R , T cSA = 6200(50) K, M B = 1.368(16) M , R B = 1.584(20) R and 
TeffB = 6200(50) K. The metal content of the system is determined to be [Fe/H] = 
+0.08(15) dex. 

In addition, the eclipsing binary has a K-type close visual companion at about 
10-arcsec separation, which is shown to be physically linked, thus sharing a common 
origin. The effective temperature of the visual companion {T e sc = 5250(200) K) is 
determined from synthetic spectrum fitting, and its luminosity (logi/L = —0.27(6)), 
and therefore its radius (R = 0.89(9) R ), are obtained from comparison with the 
apparent magnitude of the eclipsing pair. 

The observed fundamental properties of the eclipsing components are compared 
with the predictions of evolutionary models, and we obtain good agreement for an age 
of 2.6 Gyr and a chemical composition of Z = 0.026 and Y = 0.26. Furthermore, we 
test the evolutionary models for solar-mass stars and we conclude that the physical 
properties of the visual companion are very accurately described by the same isochrone 
that fits the more massive components. 

Key words: stars: individual: CD Tau - binaries: eclipsing - binaries: visual - stars: 
evolution - stars: abundances - stars: late-type. 



1 INTRODUCTION 



Double-lined eclipsing binaries are the only objects, apart 
from the Sun, for which fundamental and simultaneous de- 
terminations of masses and radii can be obtained. These 
determinations are possible through the analysis of spectro- 
scopic data in the form of radial velocity curves, and from 
modelling the photometric data in light curves. In addition, 
if we consider detached double-lined eclipsing binaries, no 
significant mass transfer has occurred between the compo- 
nents, since they are smaller than their respective Roche 
lobes. In such a case, the mutual interaction can be safely 
neglected and the components can be assumed to evolve 
like single, individual stars. Therefore, these objects yield 
simultaneous absolute dimensions for two single stars that 
are supposed to have a common origin both in time and 
chemical composition. In this situation, evolutionary models 



should be able to predict the same age for both components 
for a certain chemical composition. 

Eclipsing binaries that are members of physically bound 
multiple systems provide additional constraints for the anal- 
ysis of evolutionary models. All the information that can be 
extracted for the companion/s (effective temperature, mag- 
nitude difference, mass, . . . ) should be also fitted by the 
same isochrone that fits the eclipsing binary pair. Not many 
studies have made use of this possibility up to date, mainly 
because of the scarce amount of information available for 
the additional stellar companions. 

CD Tau (HD 34335, HIP 24663) is a bright eclipsing 
binary (Knax = 6.75) composed of two similar F6 V stars, 
which was suggested by Chambliss (1992) as a possible mem- 
ber of a triple system. It has a K-type close visual companion 
(CD Tau C) at 9.98 arcsec (ESA 1997) with a magnitude 
of V = 9.9. The angular distance between the eclipsing pair 
(CD Tau AB) and CD Tau C is large enough to allow the 
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possibility of obtaining relevant astrophysical information 
for all three components separately. 

Evolutionary models for stars with M ^ 1.2 Mq are 
known to provide a reasonably good description of the ob- 
served stellar properties, at least within the main sequence. 
The relatively large amount of eclipsing binary and open 
cluster data in this mass range (Andersen 1991; Pols et 
al. 1998) has been of remarkable help for reaching such 
an agreement. Nevertheless, recent studies (Popper 1997; 
Clausen, Helt & Olsen 1999) have pointed out that current 
stellar evolutionary models are unable to appropriately de- 
scribe the observed properties of low-mass eclipsing binaries 
(0.7-1.1 Mq). The authors have found that the models for 
the less massive stars seem to predict radii that are too 
small (log g too large) and effective temperatures that are 
too high when compared to observational data. Clausen et 
al. (1999) suggest several possible causes of the discrepancy: 
observational problems, lack of chemical composition deter- 
minations and the inaccuracies of the physics of the models, 
in particular the mixing length parameter. In any case, the 
number of systems in this mass range is still very small to 
reach statistically significant conclusions. CD Tau C, which 
is apparently a sub-solar mass star, may yield useful infor- 
mation for evolutionary model testing and could also help 
in the solution of the dilemma indicated by Popper (1997) 
and Clausen et al. (1999). 

In this paper, we give a detailed description of the avail- 
able observational data on CD Tau and their analysis for the 
determination of the fundamental properties of the three 
components. Moreover, the existence of a physical link be- 
tween CD Tau AB and CD Tau C is proved. We finally com- 
pare the resulting parameters with the evolutionary model 
predictions in order to assess their ability to reproduce the 
observations. 



2 MODELLING THE RADIAL VELOCITY 
AND LIGHT CURVES 

A radial velocity curve of CD Tau was published by Pop- 
per (1971), but the quality of the data and the poor phase 
coverage do not allow an accurate determination of the 
masses. Therefore, new radial velocity observations were 
secured with the CORAVEL scanner (Baranne, Mayor & 
Poncet 1979) attached to the 1-m telescope of the Obser- 
vatoire de Geneve and located in the Observatoire d'Haute 
Provence (France). Two different observational campaigns 
(1997 February and October) allowed us to collect 48 and 46 
measurements for the primary and secondary components, 
respectively, and 4 measurements for the visual companion. 
Heliocentric corrections to both Julian date and radial ve- 
locity were applied, and the correlation dips were analysed 
by considering colour dependence (see Baranne et al. 1979; 
Duquennoy, Mayor & Halbwachs 1991). The final radial ve- 
locity measurements for all components, together with their 
errors, are listed in Table |l| 

The solution of a radial velocity curve for a circular 
orbit, such as that of CD Tau AB, is a conceptually very 
simple problem. The only free parameters are the two ve- 
locity semi-amplitudes (Ka and Kb), the systemic veloc- 
ity (7) and a phase offset. The ephemeris for CD Tau, for 
both radial velocity and light curve solutions, was adopted 



Table 1. CORAVEL radial velocity observations of CD Tau A, 
B and C. 
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Figure 1. Radial velocity curve fit to CD Tau AB CORAVEL 
observations from 1997 February and October. Circular orbit is 
adopted. The r.m.s. residual of the fit is 0.3 km s _1 . 



as Mini = 2441619.4075 + 3? 435137 E (Kholopov et al. 
1987). We used the SBOP program (created by Dr. P.B. 
Etzel in 1978 and later revised several times) and adopted 
the Lehmann-Filhes method (Lehmann-Filhes 1894; Under- 
hill 1966) for simultaneous solution of a double-lined radial 
velocity curve. The parameters resulting from this analysis 
are shown in Table g and a plot of the observed radial ve- 
locities and the best-fitting curve is presented in Fig. [lj. The 
r.m.s. residual of the fit is as small as 0.3 km s _1 . The errors 
of the parameters were conservatively adopted as twice the 
standard errors provided by the SBOP program. Our results 
show good agreement with the analysis of the previously 
published radial velocity curve of Popper (1971), although 
the curve coverage and the individual accuracy of the mea- 
surements are significantly better in our study, leading to 
smaller formal errors. 

Three complete photoelectric light curves of CD Tau 
have been observed to our knowledge by Srivastava (1976), 
Wood (1976) and Giilmen et al. (1980). The best data are 
those of Wood, obtained in the Stromgren system, but, un- 
fortunately, the individual measurements were not published 
and are not available. Thus, for this work, we made use of 
the observations of Giilmen et al. with 292 measurements in 
each of B and V filters. 

The analysis of light curves was done by using an im- 
proved version of the Wilson & Devinney (1971, hereafter 
WD) program (Milone, Stagg & Kurucz 1992; Milone et 
al. 1994) that includes Kurucz (1991, 1994) ATLAS9 atmo- 
sphere models. As suggested by the uncomplicated shape of 
the light curve (Fig. El), we chose a detached configuration 
with coupling between luminosity and temperature when 
running the solutions. Reflection and proximity effects were 
considered for the sake of completeness, but their effect is 
expected to be negligible for such a well-detached system. 
The bolometric albedos were set to 0.5 (typical for convec- 
tive atmospheres) and the gravity darkening exponents were 
adopted as 0.28 and 0.30 (A and B component, repectively) 
as deduced from the new computations of Claret (1998). The 
mass ratio q was fixed to the spectroscopic value of 0.948, 



circular orbit was adopted and the temperature of the pri- 
mary star was set to 6200 K (see Sect. ^). Finally, both 
components were assumed to rotate synchronously with the 
orbital period. 

Both BV light curves were solved for simultaneously in 
order to achieve a single, mutually consistent solution. The 
free parameters in the light curve fitting were: the phase off- 
set, <f>o, the inclination, i, the temperature of the secondary, 
Tb, the gravitational potentials, Qa and J1b, the monochro- 
matic luminosity of the primary, La, and a third light con- 
tribution, F3. 

An automatic procedure was employed to run the WD 
program. At each iteration, the differential corrections were 
applied to the input parameters to build the new set of 
parameters for the next iteration. The result was carefully 
checked in order to avoid possible unphysical situations (like 
e.g. Roche lobe filling in detached configuration). A solution 
was defined as the set of parameters for which the differen- 
tial corrections suggested by the WD program were smaller 
than the standard errors during three consecutive iterations. 
However, when a solution was found, the program did not 
stop. Instead, it was kept running in order to test the sta- 
bility of the solutions, to evaluate their scatter and to check 
for possible spurious solutions. In the case of CD Tau, sev- 
eral runs done from different initial conditions converged to 
a stable minimum. We finally adopted the parameters that 
are listed in Table ^. 

A especially conflictive parameter is the third light since 
possible strong correlations with other free parameters may 
drive us to a wrong solution. The correlation matrix was 
checked for this possibility and we noticed weak correlations 
between F3 and all the remaining parameters but the incli- 
nation, for which the correlation coefficient was close to 0.9. 
This is a predictable result since the presence of third light 
affects the relative depths of the eclipses. Actually, the in- 
clusion of F3 as a free parameter in the solution of the light 
curve should only obey to an external evidence of some light 
excess. For this system, the situation is clear: the photomet- 
ric observations of Giilmen et al. (1980) included the vi- 
sual companion inside the diaphragm so that some amount 
of extra light is certainly present. Since we know that the 
V magnitude difference is about 3 mag and that the spec- 
tral type is later than that of the eclipsing pair components 
(A(B — V) ~ 0.4, from Tycho photometry in the Hippar- 
cos and Tycho catalogues, ESA 1997), we expect an excess 
of light of about 6.3 per cent in V and 4.4 per cent in B. 
These figures are actually very close to the values derived 
from the light curve analysis, and the colour behaviour is 
also reproduced, thus giving additional physical meaning to 
the best-fitting parameters obtained. 

The uncertainties in the parameters presented in Table 
|2] were carefully evaluated by means of two different ap- 
proaches. On the one hand, the WD program yields the 
standard error associated to each adjusted parameter. On 
the other hand, the error can also be estimated as the r.m.s. 
scatter of all the parameter sets corresponding to the it- 
erations between the first solution and the last iteration 
(about 100). We shall mention that this is, in general, larger 
than the r.m.s. scatter of the solutions alone, i.e. those that 
fulfill the criterion previously described. We conservatively 
adopted the uncertainty as twice whichever of the two error 
estimates was the largest. The r.m.s. scatters of the residu- 
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Table 2. Orbital and stellar parameters for CD Tau A and B derived from the analysis of the light and radial velocity curves, standard 
photometry and spectroscopy. The numbers in parentheses are the adopted errors affecting the last digits of each parameter. 



Radial velocity curve solution 



K A = 96.8(5) km s" 1 «j = |ia = 0.948(7) 
K B = 102.1(5) km s" 1 a = 13.52(7) R 



7 = -29.3(3) km s" 1 
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F 3 \ B = 0.052(15) b 
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Photometry and spectroscopy 



T eS = 6200(50) K 
frotA = 28(3) km s _1 



[Fe/H] = +0.08(15) dex 
frotB = 26(3) km s _1 



E(B - V) = 0.0 mag 



Stellar physical properties 



M A = 1.442(16) M 


M B = 1.368(16) M Q 


R A = 1.798(17) Rq 


R B = 1.584(20) R 


log pa = 4.087(10) dex 


logg B = 4.174(12) dex 


T effA = 6200(50) K 


T effB = 6200(50) K 



Adopted, b Phase dependent quantity: reference phase = 0P25 
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Figure 

2. Light curve fit to the observed B and V Comparison— Variable 
(Cp— V) differential photometry of CD Tau. Also shown are the 
Observed— Computed (O— C) residuals. Solutions obtained with 
WD program applied to the observations of Giilmen et al. (1980). 



als in the light curve fit are 0.011 mag (n — 292) for both B 
and V. Fig. g shows the light curve fit to the observed B and 
V differential photometry and the corresponding residuals, 
where no systematic trends are observed. 



We shall compare our analysis with that of Wood 
(1976), Giilmen et al. (1980) and Russo et al. (1981). Re- 
garding the first two works, significant differences exist in 
the inclination and in the fractional radii. The inclination is 
0.8-1° larger in our solution, r\ is also 0.010 larger in our 
case, but re is about 0.005 smaller, the latter implying a 
large change in the ratio of radii. However, our computed 
elements are in very good agreement with those of Russo 
et al. (1981). They also employed the WD program (al- 
though an earlier version) to analyse the light curves, while 
the other authors used Wood's model. Thus, the reason for 
the discrepancy in the solutions is probably related to in- 
trinsic differences of the models. The radii of Wood (1976) 
and Giilmen et al. (1980) lead to a physically unlikely situa- 
tion for a detached system, with a secondary component too 
evolved when compared to the primary (log qa <; log <;b). 

As a brief summary, our final solution shows a system 
composed of two similar stars of about 1.4 Mq in a moder- 
ately evolved stage (log g ss 4.1), although still in the main 
sequence. The remarkable point is that the absolute dimen- 
sions of the components are known to a precision of about 
1 per cent. 



3 EFFECTIVE TEMPERATURE 
DETERMINATION 

Stromgren-band and IR-band photometry were used for 
the determination of the "mean" effective temperature of 
CD Tau AB. Standard Stromgren indices for CD Tau AB 
were taken from the catalog of Hauck & Mermilliod (1998) 
and are shown in Table ^. We adopted the extinction cor- 
rection calibration of Crawford (1975) and the photometric 
grids of Napiwotzki (1998), based on Kurucz ATLAS9 at- 
mosphere models. The application of these calibrations in- 



© 1999 RAS, MNRAS 000,|l|-| 



CD Tau: a detached eclipsing binary with a solar-mass companion 5 



Table 3. Stromgren (from Hauck & Mermilliod 1998) and IR-band (this work) observations for CD Tau AB. The JHK measurements 
were taken at quadrature, when full light contribution from both components is present, n is the number of averaged photometric 
measurements. 
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dicated a negligible amount of interstellar extinction and 
yielded a value of 6200 K for the effective temperature, as 
shown in Table |§| 

In order to better constrain the effective temperature 
determination, several measurements of CD Tau in the IR 
bands were obtained in the 1 m Carlos Sanchez Telescope 
(Tenerife, Spain). The observations were made on 1998 
February 19 and 24, and April 1 and 3, when the eclips- 
ing component was near quadrature. The reduction of the 
data and the transformation to the J, H and K standard 
filters was done by using the procedure explained in Masana 
(1999). The Infra-Red Flux Method (IRFM, Blackwell et al. 
1990) was employed to obtain an accurate determination of 
the effective temperature, by comparing the observed flux 
distribution with Kurucz ATLAS9 model atmospheres. A 
log g of 4.1 and [Fe/H] = +0.08 (see Sect. 0) were adopted 
as input parameters in the analysis, although their influence 
is almost negligible. The resulting values of the effective tem- 
perature obtained at each passband are listed in Table |^. 

As it can be seen, the mutual agreement between the 
three JHK determinations is very good, and they are almost 
identical to the independent value obtained from Stromgren 
photometry. Therefore, we can confidently adopt a T e g = 
6200±50 K, where the error includes both the variance of the 
average for all filters and the calibration error of the IRFM 
method (see Alonso, Arribas & Martinez- Roger 1996). 



4 CHEMICAL COMPOSITION 

Two high-resolution spectra (0.15 and 0.12 A/pix) in differ- 
ent spectral regions were obtained with the aim of determin- 
ing the atmospheric chemical composition of CD Tau AB. 
The instrument used was an echelle spectrograph attached 
to the 74-inch (1.9 m) telescope of Mount Stromlo Observa- 
tory (Australia). The spectra were obtained on 1997 Octo- 
ber 17 and 19, when the eclipsing component was at phase 
0?87 and 0?44, respectively. One of the spectra, with an in- 
tegration time of 900 s, spans 5730 to 6045 A (around the 
Nald line of the interstellar medium), whereas the other one 
covers the range from 3830 to 4070 A (including the strong 
lines of Can K and H) with an integration time of 600 s. 
The S/N ratio is of the order of 60 and 20, respectively. As 
a preliminary reduction, the raw images were corrected for 
bias, dark current and flat field using the standard proce- 
dures (see Jordi et al. 1995) and then wavelength calibrated. 
Radial velocity corrections were applied to refer the wave- 



length scale to the heliocentric reference frame and, finally, 
the flux scale was normalized to unity. 

The element abundances were determined by fitting of 
synthetic spectra to our observed spectra. Unfortunately, 
their S/N ratios (especially for the one around Can lines) 
are too low to allow a very accurate determination of the 
chemical composition (see Cayrel de Strobel 1985; Edvards- 
son et al. 1993) but still enough for a preliminary estimation. 

For the synthetic spectrum fitting, we adopted Kurucz 
ATLAS9 models and the set of programs developed by Dr. 
I. Hubeny (synspec, rotin and synplot). Nevertheless, 
we modified these programs by considering the double-lined 
nature of our spectra, in such a way that they use two at- 
mosphere models, two rotational velocities, two turbulent 
velocities and two radial velocities as input data. The two 
resulting synthetic spectra are subsequently co-added. Since 
we are dealing with an eclipsing binary system, the sur- 
face gravities and the temperature ratio obtained by mod- 
elling the light and radial velocity curves provide us with 
useful constraints. The "mean" system log g is 4.1 dex, so 
we adopted Kurucz atmosphere models at a log g of 4.0 dex. 
Nevertheless, several tests with log g — 4.5 dex showed a 
very weak influence of this parameter on the synthetic spec- 
tra. The effective temperatures are strongly correlated with 
the atmospheric chemical composition determination. Thus, 
it is necessary to adopt accurate and mutually consistent 
effective temperatures. From the very precise T e g determi- 
nations through IR and Stromgren photometry, we adopted 
atmosphere models computed at T e g = 6200 K for both 
components. Since the temperatures are identical, but not 
the luminosities (i.e. the radii), a factor of 0.77 (see Table ^ 
was applied to the flux output of the secondary component 
when co-adding the synthetic spectra. 

The radial velocity shifts of each component were 
computed from the radial velocity curve at the ac- 
tual phases and introduced as a fixed parameter 
(42 km s _1 and —103 km s -1 for the Na id-region spectrum 
and —65 km s" 1 and 8 km s _1 for the Call-region spec- 
trum). The rotational velocities of each component were ini- 
tially adopted as the synchronization values, 26 km s _1 and 
23 km s" 1 , and were subsequently changed to 28 km s _1 and 
26 km s _1 in order to obtain the best fit to the line profiles. 
The error of these rotational velocities is estimated to be 
about 3 km s _1 . 

Taking into account all the adopted input parameters, 
the spectrum was studied in small intervals 20 A wide. From 
several fitting tests at different Fe-peak element abundances 
(-0.1 dex, -0.05 dex, 0.0 dex, +0.05 dex, +0.1 dex and 
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Figure 3. Fit to an echelle spectrum of CD Tau AB (double- 
lined, phase = 0?44) from 3910 to 4070 A obtained at Mt. Stromlo 
(Australia) on 1997 October 19. The thick line represents the ob- 
served spectrum whereas the thin line is the best-fitting synthetic 
spectrum computed using Kurucz ATLAS9 fluxes. The parame- 
ters used for the synthetic spectrum computation as well as some 
comments to individual features can be found in the text. 

+0.15 dex) we deduced that the observed spectrum is best 
reproduced with a slightly super-solar composition. We esti- 
mate a value of [Fe/H] = +0.08±0.15 dex. The lack of lines 
of other elements does not allow abundance determinations, 
but there exists evidence of a slight underabundance of Na 
and a slight overabundance of Si, however, with low reliabil- 
ity. The parameters derived from the spectrum analysis are 
listed in Table ^. 

In general, the spectrum covering 5730-6045 A was 
found to be of little help to obtain a Fe abundance determi- 
nation since there are no strong lines of Fe-peak elements, 
although we obtained a good fit to the stellar Nald line. 
Moreover, the reddest zone of the spectrum (A > 5970 A) 
was strongly contaminated by non-stellar features. 

The spectrum spanning 3830 to 4070 A (which is 
mostly shown in Fig. []) proved to be much more use- 
ful for abundance determination. The wavelength region of 



A < 3870 A showed poor flux detection statistics and it was 
not used in the fit. The group of Fe, Mg and V lines at 
A3878 A were fitted reasonably well, although the blue edge 
has a slightly different slope. The Balmer line at A3889 A is 
somewhat deeper in the synthetic spectrum, but this may 
be caused by a not very detailed treatment of this broad line 
in the latter. The lines of Fe-peak elements between A3890 
and A3910 A can be fitted perfectly. However, although the 
depth of the Fe lines in the 3918-3920 A region is well re- 
produced, there seems to be a problem with the continuum 
level, which shows a dip in the observed spectrum. 

The fit to the Can K (A3933.7 A) and H (A3968.5 A) 
line profiles is very good and was useful for determining the 
rotational velocities along with some strong Fe lines. From 
the Ca n H line to A4040 A the fit to the observed absorption 
lines is very satisfactory in such a way that the model and 
the observations are hardly distinguishable in most parts. It 
is only between A4020 and A4028 A that a slight underab- 
sorption is present in the observed spectrum, but it seems to 
be related to some continuum effect. In the remaining zone 
of the spectrum, we will only point out a couple of features 
besides the general good agreement. The Fel line profile at 
A4046 A presents a subestimated absorption, especially for 
the primary component. This behaviour, not observed in any 
other strong line, appears to be caused by the contribution 
of another absorption source. This is the strongest Fel line 
in this region, and some amount of interstellar absorption 
at v r ~ may be present so that it would be superimposed 
on the primary component line. The check of such a hy- 
pothesis would require the analysis of some spectra of the 
same region at different phases. Finally, the absorption Fel 
feature around A4055 A is overestimated in the synthetic 
spectrum while another quite strong Fe I line at A4064 A is 
well reproduced. 

In conclusion, the agreement between the observed and 
the synthetic spectrum computed with the parameters listed 
in Table ^ is good, as shown in Fig. ^[ In some features, the 
slight flux underestimates may be caused by the poor S/N 
of the observed spectra, by some deficiencies in the synthetic 
spectra or, more unlikely, in the model atmosphere fluxes. 



5 THE VISUAL COMPANION: CD TAU C 

Radial velocity measurements of CD Tau C were also ob- 
tained with the CORAVEL scanner. In spite of the faintness 
of this star, near the limiting magnitude of the instrument, a 
very clear correlation dip was present due to the closeness of 
its spectral type to that of the mask (K2 III). The four mea- 
surements lead to a mean v T = —27.8 ± 0.8 km s _1 , which 
should be compared to the systemic velocity of the eclipsing 
system, 7 = —29.3 ± 0.3 km s~ , in order to assess the pos- 
sibility of a physical link. As it is seen, the agreement is fair. 
Moreover, the magnitude difference between CD Tau AB 
and CD Tau C is compatible with what is expected for a 
K-type star located at the same distance as the eclipsing 
pair. 

From Hipparcos measurements, the distance to CD Tau 
is known to be 73±9 pc, which means that an angular 
separation of 9.98 arcsec can be translated into 730 AU. 
Assuming the total mass of the system to be around 
3.7 M Q (0.9 M Q is adopted for CD Tau C) and the current 
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distance to be the semi-major orbital axis, Kepler's Third 
Law leads us to an orbital period of fa 10 000 yr and a max- 
imum expected radial velocity difference of « 2.1 km s _1 . 
Such a wide system seems not to be capable of surviving 
any weak encounter with another object. Thus, although it 
cannot be assured that CD Tau AB and CD Tau C are in or- 
bit around each other, there exist evidences suggesting that 
they actually form a common pair. 

The determination of the largest possible number of pa- 
rameters for the visual component may be extremely use- 
ful when evaluating the model predictions. Unfortunately, 
neither the precision of the Stromgren photometric indices 
nor the intrinsic accuracy of the calibrations is enough 
for a reliable determination of the effective temperature of 
CD Tau C. Moreover, no IR observations could be collected 
due to the strong contamination caused by the bright com- 
panion. Thus, no reliable temperature estimation was pos- 
sible from a photometric point of view. 

However, an spectroscopic temperature determination 
of CD Tau C was performed via synthetic spectrum fitting. 
Along with the spectra of CD Tau AB, two spectra of its 
faint visual companion were also secured, with identical in- 
strument configurations. The S/N ratios of the spectra were 
of the order of 4 (the star is about 3 mag fainter), so that 
any possibility of abundance determination was ruled out. 
However, since there are strong evidences that CD Tau AB 
and CD Tau C are a physical system, we can reasonably as- 
sume that they share the same chemical composition. Thus, 
the spectra can provide a valuable estimation of the effective 
temperature of this early K-type star, for which there is no 
other clue. 

The spectrum fitting was performed in a similar way as 
for CD Tau AB, but taking into account the single nature 
of the star. The surface gravity was adopted as 4.5 dex, the 
radial velocity as —28 km s^ 1 (from CORAVEL observa- 
tions) and the same chemical composition as CD Tau AB. 
Synthetic spectra were computed from Kurucz models at 
250 K intervals between 4750 K and 5500 K and compared 
to the observed spectrum. In many zones, the choice of the 
best-fitting temperature was not clear due to the noise, but 
there were two regions (around Ca II H line, which shows 
some chromospheric emission, and around A4045 A, with a 
strong Fe I line) that allowed us to adopt 5250 K as the 
effective temperature that yields the best fit. The error of 
this determination is estimated to be around 200 K. 

Further information can be extracted from the mag- 
nitude difference between the eclipsing pair and the visual 
component. In the Hipparcos catalogue, the following val- 
ues of the Hp magnitude are listed for CD Tau AB and 
CD Tau C: 

Hpab = 6.871 ± 0.007 Hp c = 9.932 ±0.110 

that lead to a magnitude difference AHp — Hpc — Hpab = 
3.06 ± 0.12 mag. From AHp we derive: 

AMboi = AHp + ABC(Hp) 

where BC(Hp) is the bolometric correction for the Hp pass- 
band. Straightforward transformations lead to: 

From the BC calibration of Cayrel et al. (1997) we adopted 
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Figure 4. H-R diagram of stars with masses comprised between 
0.75 and 1.25 Mq. Empty symbols are primary and secondary 
components of eclipsing binaries, whereas the filled symbol rep- 
resents CD Tau C. 

BC(H p )ab = -0.15 ± 0.05 and BC(Hp) c = -0.35 ± 
0.05, leading to a ABC{Hp) = -0.20 ± 0.07. These val- 
ues and their estimated errors can be combined to derive 
Alog(L/L ) = -1.14 ±0.06. 

The absolute luminosities of CD Tau A and B can be 
computed from their radii and effective temperatures and 
we obtain: 

log(L/L Q ) A = 0.63 ± 0.02 log(L/L ) B = 0.52 ± 0.02 

that can be combined to compute the total luminosity of 
the eclipsing pair, which comes out to be ^(Lab/Lq) = 
0.87 ± 0.02. We finally obtain the luminosity of CD Tau C 
that is: 

log(L c /L Q ) = -0.27 ±0.06 

Also, from the effective temperature and the luminosity, 
the radius can be easily computed and we derive Rc = 
0.89 ± 0.09 Kq. All these values are also in agreement with 
the results obtained when adopting the magnitude difference 
listed in Wood (1976). 

For illustrative purposes, we present in Fig. ^| an H- 
R diagram showing the position of CD Tau C and a sam- 
ple of well-studied eclipsing binaries in the mass range 
0.75 < (M/M Q ) < 1.25, taken from the compilation of 
Clausen et al. (1999). As it can be seen, CD Tau C is located 
within the distribution described by all the eclipsing binary 
components. 



6 COMPARISON WITH EVOLUTIONARY 
MODEL PREDICTIONS 

The comparison of model predictions and the observed prop- 
erties of CD Tau AB was done under the assumption that 
the evolutionary models should be able to fit an isochrone 
to both components of the system for a certain chemical 
composition. We made use of the evolutionary models of 
CG (Claret 1995, 1997; Claret & Gimenez 1995, 1998) that 
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Figure 5. log g — log T e g representation of evolutionary tracks 
and isochrone for CD Tau A and B. They are computed for the 
best-fitting value of the chemical composition by using the evo- 
lutionary models of CG (see text for references). 



Figure 6. H-R diagram showing the position of CD Tau C and 
evolutionary tracks for initial masses from 0.8 to 1.0 Mq. A 
2.5 Gyr isochrone (the age of CD Tau C) is also shown in the plot. 
The chemical composition was fixed at Z = 0.026 and Y = 0.26. 



constitute a grid with several (Z, Y) values. The actual 
isochrone fitting was done by means of an interpolation rou- 
tine that uses the masses, radii and effective temperature of 
the components and yields the best-fitting values for the age 
and the chemical composition (Z, Y) of the system (among 
other parameters). Thus, both Z (initial metallicity) and Y 
(initial helium abundance) were treated as free parameters. 
The computational details of this algorithm will be provided 
in a forthcoming paper (Ribas et al. 1999). 

The evolutionary models predict an age of 2.6 ±0.3 Gyr 
and best- fitting chemical composition values of Z — 0.026 ± 
0.007 and Y = 0.26 ± 0.04. Indeed, the metal abundance 
derived from the isochrone fitting is in good agreement 
with the spectroscopic metallicity determination: [Fe/H] — 
+0.08T0.15 or Z = 0.023T0.008. A log ff-log T cff represen- 
tation of CD Tau A and B, together with the evolutionary 
tracks and isochrone computed for the best-fitting chemical 
composition, is presented in Fig. 

Considering the physical relation between CD Tau AB 
and CD Tau C, a common origin, both in time and composi- 
tion, can be assumed. This implies that the same isochrone 
should simultaneously fit all three stars at an age of about 
2.6 Gyr. Due to the presumably large mass difference be- 
tween the components, it is advisable to correct this esti- 
mated age by taking into account the different lifetime of 
the stars in the pre-main sequence (PMS) phase. The evo- 
lutionary tracks published by D'Antona & Mazzitelli (1994) 
indicate that the difference in duration of the PMS phase 
(the mass of CD Tau C is supposed to be 0.8-0.9 Mq) is 
about 0.1 Gyr. So, the age of CD Tau C has to be corrected 
to a value of 2.5 Gyr when comparing to post-ZAMS evolu- 
tionary models. 

In order to compare the observed data with the model 
predictions, evolutionary tracks of initial masses 0.8, 0.85, 
0.9, 0.95 and TO Mq for the adopted chemical composition 
were kindly provided by Dr. A. Claret. They were computed 
with the same input physics of the CG evolutionary mod- 
els. A log T c ff — log L plot of the evolutionary tracks and 



the observed values for CD Tau C is presented in Fig. |6[ 
A 2.5 Gyr isochrone is also plotted. As it can be seen, the 
theoretical predictions and the observational data show an 
excellent agreement, well within the error bars. Notice that 
no free parameters remain, since the chemical composition 
of the star is also known. The best-fitting parameters that 
the models predict for an age of 2.5 Gyr are a mass of about 
M = 0.99 M and a log g of 4.54 dex. 

Further refinement of the effective temperature of 
CD Tau C (by means of e.g. high-accuracy IR measurements 
of the visual component alone) may provide a more strin- 
gent test to the quality of the evolutionary models. What 
seems clear at this point, is that the systematics that Pop- 
per (1997) and Clausen et al. (1999) pointed out are not 
observed in CD Tau C. Indeed, the stellar parameters are 
very accurately described by the isochrone that fits the more 
massive components. 
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